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ABSTRACT 

The purpose of this work is to describe a new method for calculating the aerodynamic performance of 
horizontal axis wind turbines. The method, entitled the helical vortex method, directly calculates the 
local induced velocity due to helical vortices that originate at the rotor blade. Furthermore, the 
method does not require a specified circulation distribution. 

Results of the method are compared to similiar results ol)tained from Wilson PROP code methods 
(Prandtl. Goldstein, NASA and no tip loss) as well as to existing experimental data taken from the NASA 
Mod-0 wind turbine. It is shown that results of the proposed method agree well with experimental values 
of the power output both near cut-in and at rated wind speeds. Further, it is found that the method does 
not experience some of the numerical difficulties encountered by the PROP code when run at low wind 
velocit ies , 


INTRODUCTION 

Recently there has been a renewed interest 
in wind turbines as a means of producing power. 
Considerable effort has been directed toward 
exper imental and theoretical studies of wind 
turbine performance. 

At present, the analytical method used in 
predicting aerod^^namic wind turbine performance by 
the NASA Wind Energy research group utilizes _.a 
modified blade element theory [1] or Glauert vortex 
thecFry incorporated with Prandtl [2] , Goldstein 
[3] and NASA [4] tip loss models. In the Glauert 
vortex theory, it is assumed that trailing vortices 
originate from the rotating blades and form a 
helical vortex system that passes down stream. 

This vortex system, in turn, induces velocities 
which alter the flow around the blade. However, 

In this method, the induced velocities are not 
deteimilned directly due to the complexity of the 
method. instead, an interference velocity is 
calculated as the induced velocity of this vortex 
system on any blade clement. The calculation of 
these interference velocities is simplified by 
assuming that the rotor has an infinite number of 
blades. The latter assumption removes the com- 
plexity associated with the periodicity of the 
flow and permits mom ent um theory to be directly 
used to evaluate the interference velocities. 
However, for single and double bladed rotors, this 
assujaption may be inadequate. 

Consider a wind turbine, rotating with an 
angular velocity about its horizontal axis of 
rotation, that is placed in a uniform stream of 
wind of velocity parallel to the axis of 
rotatiori. The blaae can be coned or tilted to an 
angle of y from the plane perpendicular to the 
axis of rotation. The velocity experienced by the 
typical blade is shown in Fig. 1. It can be seen 
from Fig. 1, that Vo(l-a)cosi^ is the wind free 
stream velocity minus the axial interference 
velocity normal to the blade surface. Further, 

Qr fl +a * )cos?|; is the relative velocity of blade to 
the air velocity and accordingly is the angular 
velocity of the blade section plus the rotational 


interference velocity. The factors, ^ and ^ are 
the axial and the rotational interference factors 
respectively. 


In order to evaluate the local drag force 
and torque, the interference factors _a and ^ must 
first be determined. However, before presenting 
the equations for determining ^ and a_|_, it may be 
helpful to review the tip and hub losses for a 
blade. It should be anticipated that the flow 
around a rotor blade of finite length will be 
disturbed at the tip and root or hub of the blade. 
These disturbances are due to the fact that the 
difference of pressure between the upper and lower 
sides of the blade disappear along the blade span 
ond must therefore diminish to zero toward the tip 
and hub. An approximate method of estimating the 
effect of these radial disturbances has been given 
by Prandtl [2]. Later, a more accurate solution 
to this problem was developed by Goldstein [3] . 
Introducing a reduction factor, F, to account for 
the fact that only a fraction of the air between 
successive vortex sheets of the slip stream 
receives, the full effect o f the motion of these 
sheets, into the momentum equation for the flow 
at radius r and then equating the differential 
drag and torque equations obtained respectively 
from momentum and blade element theories, permits 
the writing of equations which must be satisfied 
by a and a^ , These equations may be written as; 


a o^iC^cos^ + Cj^sin(j>)cos^\|^ 
1-a 4 F 


( 1 ) 


and 

a- 

1+a' " 4 F sinsj) cos<j> 


( 2 ) 


Expressions for the reduction factor F can be 
found in [2,3]. The tip loss method which has 
been used by NASA is the so-called ’’effective 
radius model”. Due to the radial flow near the 
boundary of the slipstream, there is a drop of 
circulation which can be represented by an equiv- 
alent rotor with an infinite number of blades 
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and with the same drag force but with a smaller 
radius R^, i.e,, an effective radius is defined as 


where B is a constant tip loss factor, 
o ^ 

Wilson and Walker [4] developed a computer 
program called PROP to evaluate the solution of 
eqs. (1) and (2). A preassigned tip loss factor 
(F = 1, no tip loss) was used and an iteration 
procedure for determining a and a' for a given 
differential element was employed. Once a and a’ 
were determined, torque, drag and power could aTT 
be calculated from the appropriate equations. 
However, the value of the axial interference factor 
a has a limit. For a positive axial velocity V , 
the induced axial velocity of out flow near the^ 
rotor with a coning angle is (V a cos;//), thus, 
from vortex theory, the induced axial velocity in 
the ultimate slipstream will be (2V a cosi//) . 
Therefore, the trailing velocity of air far behind 
the rotor is [V^(l-2a)cos\{/] . Equations (1) and (2) 
are based on the notion that the axial velocity at 
any location is unidirectional, and, therefore, the 
analysis is valid only if a £ J , If the rotor 
absorbs all the energy, i.e., = 0, then a would 

have a maximum value of 1 . Unfortunately, by 
using the PROP code to calculate a, under certain 
operating conditions, particularly for large tip 
speed ratios (small wind velocities), the local 
value of the axial interference factor can exceed 
this limiting value for all tip loss models. 
Therefore, the equations cease to be valid. It 
was also found that the PROP code generally under- 
estimated the performance (5, 6] , Accordingly, it 
is the purpose of this paper to develop a computa- 
tional method for predicting aerodynamic perfor- 
mance of a horizontal axis wind turbine that avoids 
these difficulties. In this method, the induced 
velocity is directly calculated by integration of 
Biot-Savart ' s law under the assumption that a 
filament of the trailing vortices is helical in 
nature, extends infinitely downstream of the rotor, 
and has constant pitch and diameter. Of primary 
importance, the method does not use interference 
factors. 

The predicted performance of wind turbines 
obtained by use of the present technique will be 
compared with those obtained from the PROP code and 
with some experimental data obtained from the Mod-0 
100 kw wind turbine operated by NASA - Lewis 
Research Center in Sandusky, Ohio. 

ANALYSIS 

General Assumptions and Discription of the Probl em 

Two major assumptions upon which the present 
analysis is based may be stated as follows: 

1. ) The trailing vortex system is helicoidal 
with constant pitch and diameter, and extends 
infinitely far downstream from the blade. 

2. ) The relative velocity of a blade element 
to the medium is identical to that in two- 
dimensional motion if reference is made to 
the relative velocity between the element and 
the medium. 


Unlike Goldstein theory, the helical vortex 
assumption does not restrict the circulation dis- 
tribution along the blade. In the Goldstein 
theory, an optimum circulation distribution was 
used that corresponded to a rigid helicoidal 
vortex system moving backward with constant 
velocity. However, in both theories, it is assumed 
that the slipstream expansion may be secondary and 
thus neglected. 

The physical model of the problem is illus- 
trated in Fig, 2. The coordinate system, which is 
shown in Fig, 3 was chosen such that the z coordi- 
nate is the distance measured from the rotor to a 
segment of the trailing vortex parallel to the 
axis of rotation of the rotor and that the r axis 
is along the blade. The conning angle of blade is 
The induced velocities are computed in terms 
of the coordinate system in Fig. 3 as will be 
demonstrated in the next section. 

Calculation of Induced Velocity 


Suppose that circulation at a particular span 
location, say r, of the blade is F. The circula- 
tion at a nearby point, say r + dr is 


r 


dr , 

dr 


The difference in these circulation values i.e.. 


dr , according to Helmholtz theory, the strength 


of a vortex that will spring out from the blade 
element of dr to form the trailing vortex, Fig. 2. 
The induced velocity normal to the resultant 
velocity W at a point r* due to the helical vortex 
filament originating from each of the blades and 
extending to infinity is determined, via the Biot- 
Savart’ law [7] by. 





C4) 


where = (^ - 4'cos0j^)4^' 

= f^(0sinOj^ + cos0^) ] A^^cos^;// 

Dj = - 255'cose^. + g^X^^cos^il, 

= X ^ cos^ij) + 

2 O r 

V 

A - ^ 

o i2Rcost// 




r 

R 




r ’ 

W 
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An expression similar to eq. (4) was employed by 
Plencner [8] for studying propeller characteristics. 

The integration with respect to 6 in eq. (4) 
is an arduous task. Furthermore, the fQrm_of the 
expression, as it stands, is inconvenient for 
practical use because the integral becomes infin- 
itely great at the point C’ = 5 and 0=0. 

Clearly, this feature unless properly handled can 
cause considerable errors when calculating the 
induced velocity. 


Moriya [7] proposed a method to remedy the 
problem. In his method, he introduced a so-called 
induction factor I which is defined as 

B - N + N 

= (5-C) I r 3/2, 1 7 2 

k=lJo ^2 


With this factor, the singularity at C 
be shov^ to reduce eq. [4) to 


K 

n 




I 


47rR cosX CC-f/) 


4 ' can 


(6) 


the value of % is usually small, it may be 
assumed that 


w 

a 

i W j 


^ d' 

l( g.g'.Ao) dg 


4ttS 7R^ cos^ip 


a - 


( 9 ) 


^hub 


The circulation function, r, may be expressed as 
a Fourier-Sine series satisfying the following 
boundary conditions 

'■(5 = O = = 0 

An expression that meets these conditions may be 
written as 

rCO = Z A^sinmw (10) 

m=l I hub 

The relation between circulation and lift coeffi- 
cient is given by the Kutta- Joukowski theorem for 
a two-dimensional airfoil as 

r = 1 c W (11) 


it should be mentioned that the induction 
far tor, which is a continuous function of X^, 

4 and 4’, is simply the ratio of the normal 
Induced velocity for a helical trailing vortex 
system to the normal induced velocity for a 
straight trailing vortex system of the same 
strength. At a point f =4, the induction factor 
I has a value of unity. In the following, the 
induction factor will be determined by numerical 
integration using 24 terms of a Lagurre-Gaussian 
Quadrature [9] , 


Goverr 


Equation Formulation 


The governing equation used in this study can 
be established from the fact that if the blade 
section is set at a geometrical angle of attack, 

IX, , relative to the incoming wind velocity, W, as 
shown in Fig. 4, its setting relative to W is 

a diminished by the downwash angle, a.. This 
g ' , . 1 

* relation may be written as 


a 


g 


a 

e 


where 



a - = tan 


w 


n 


W 


(7) 


and a is an effective angle of attack. The 
negative sign appearing in eq. occurs because 
the induced angle of attack, a., is itself a 
negative number. The geometric angle of attack 
can be obta_ined from Fig. 4 as 

a = tan. ^ B (8) 

^ 41 

where gis the pitch angle of the blade element at 
4’* 


The induced angle of attack may be expanded 
in series form in terms of Moreover, because 

W 


where c is the chord length. 


Using equations (5) and (7) through (11), the 
aerodynamic characteristics for_a given blade 
element can be determined. The calculation 
involves an iteration process and requires that 
the series for F be truncated at N terms. A 
computer program has been developed for this 
iteration and proceeds according to the following 
logic : 


1 ) As sume a 
^ e 


a in the first iteration, 
g 


2) Find from two-dimensional airfoil data 
such as that shown in Figs. (5) and (6). In 
the program, experimental data for and 

were curve fitted using polynom ial s _ 

and were stored in the program. 


3) Calculate the local circulation r(4) by 
employing eq, (11). 


4} With the values of r(4) known at each 
location 4, the constant's, A^, appearing in 
eq. (10) can be determined by solving a set 
of linear algebric equations. 


5) Calculate local ot. by using eq. (9). 

This can be done by substituting eq. (10) 
into eq. (9) and integrating the resulting 
equation. However, in practice, the integral 
is broken into three separate integrals to 
remove the singularity as follows 



1 

rV-6 

g(?)de = 

/ g(C)d? 

hub 

1 

^hub 


+ I gC4)d4 

4^-<5 



( 12 ) 


where 6 is an arbitral ly small number, say 
6 = 10'^. 
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( 16 ) 


The first and the last terms on the right hand 
side of eq. (12) can be evaluated by using stand- 
ard numerical integration procedures. The second 
integral on right hand side of eq. (12) contains 
a singularity, but, it can be shown that its value 
is very small, 0(10“'^), and, thus, may be ignored 
in the calculation. 

6) Calculate a new a - a + a.. 

e g 1 

7) Obtain from two-dimensional data using 
newly calculated 

8) Compare the current C^^ to the previous value 
of Cj^, and if equal, stop, otherwise repeat the 
iteration by returning to step 3. 

Power, Torque and Drag on Wind Turbine 

Once the iteration process is successfully 
tei'minated and the effective angle of attack 
distribution calculated, the rotor thrust, torque 
and power can be calculated respectively from the 
following equations 

^ = IpBcW^ (C^cos4i+Cj^sin<|))cosij> 

^ = ipBcW^ (C^sin<|)-CpCOS<t>)r cosi/; 


= 0.95 (Pj^ - 0.075 Pg) 
in which: 

Pg ; generated electrical power, kw 

Pj^ ; power produced by rotor, kw 

Pg ; rated electrical power, kw (100 kw) . 

Predicted alternator power is plotted against 
wind speed in Fig. 9. Also shown are correspond- 
ing theoretical results obtained by using the PROP 
computer code and existing experimental data [12]. 

From these figures, it can be seen that there 
is a closer correlation between the values pre- 
dicted by the present method and the experimental 
data than there is with the four PROP code methods 
of predictions. It is generally known that PROP 
code models underestimate performance compared with 
experimental data [5, 11, 13], However, our results 
are 10 to 15% higher than are PROP code values for 
wind speeds ranging from 8 to 20 mph. There is 
excellent agreement between our data and the exper- 
imental values. Our data predicts a "cut-in” wind 
speed of 8 mph which agrees with the actual value 
reported in [5]. It is known that the "cut-in" 
speed is not predicted accurately by any of the 
methods found in the PROP code. 


Integration yields 

T = [" (IpBcW^) (Cj cos(|)+Cj^sin4))R cos^j^dC 


hub 


Other calculations have been made for the 
NJOD-l wind turbine (2MW) . The operating conditions 
(13) and the input data for these calculations are 

listed in Table II. The airfoil used for the blade 
construction is an NACA 4418 with a "half rough" 
surface condition and a variable twist angle along 
the blade span. 


Q = / (IpBcW^) (C sin(|)-C cos 4 »)R^cos 4 ^ dC 


^hub 


(1^) The chord and twist angle are tabulated in 

Table III. 


P = / (IpBcW^) (R^Qcosij;) (C^sin4^-C^cos(J))CdC 
^hub 


(15) 


RESULTS AND DISCUSSION 


The power coefficient C vs. tip speed ratio 
for MOD-1 is presented in ^ Fig. 10. As can be 
seen, the results obtained from the PROP code are 
unlikely as the power coefficient has a saddle- 
like distribution. On the other hand, our results 
do not exhibit this behavior. It is believed that 
incorrect handling of the interference factors is 
the source of the difficulty. 


Tlie proposed method has been used to calcu- 
late performance of two-bladed wind turbines. The 
particular blade and operating conditions, listed 
in Table 1 and used in the calculations, corres- 
pond to those of NASALS Mod-0, lOOkw large wind 
turbine located in Sandusky, Ohio. Further des- 
cription of this facility may be found in [10]. 

A circulation distribution corresponding to 
the operating condition tabulated in Table 1 is 
plotted in Fig. 8. As one can see, it is quite 
irregular and departs greatly from an eliptical 
shape generally used in aerodynamic studies in 
which the distribution of circulation is unknown. 

An empirical expression, taken from [11], was 
used in the calculations to compare predicted to 
recorded alternator power. This expression 
accounts for drive train losses and is given as 
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torque coefficient 

r distance along the blade 

Tj r cos Ip 

R radius of rotor 

Re effective radius 

T thrust 

V wind velocity 

o 

w normal induced velocity 

n 

W vector sum of rotational and wind 

velocity 

W’ vector sum of W and 

z rearward distance from the rotor 

a angle of attack 

a effective angle of attack 

e 

geometric angle of attack 
induced angle of attack 
0 pitch angle 

r circulation 

p density 

BC 

o solidity = 

^ ■' L 

iji coming angle 

Q rotational speed of blade 


NOMENCLATURE 

a axial interference factor 

a' rotational interference factor 

B number of blade 

c chord 

drag coefficient 
Cj lift coefficient 

C power coefficient 

F 

L 1 ift 

P power 

0 torque 
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Table I 


Operating Condition of the Utility Pole Wind Turbine 


Root /Tip Chord 
Chord Distribution 
Percent Root Cut 
Blade Radius R, ft 
Conning Angle, deg. 

Sol idity 

Thickness to Chord Ratio 
Airfoil 
Pitch Angle 
Airfoil Surface 

Operating RPM 
Twist Angle 


6,25/2.08 
see Fig. 7 
23% 

62.5 

3.8 

0.033 

.24 

NASA 230-24 
0 

Rib Stitched, Fiberglass Cloth 
(half rough) 

32 and 33 rpm for 2 blade 
0 


Table II - Input Data and Operating Condition for MOD-1 Wind Turbine (2MW) 


Root /Tip Chord 
Blade Radius R, ft. 
Hub Radius 
Percent Root Cut 
Conning Angle, deg. 
No. of Blade 
Airfoil 
Pitch Angle 
Operating RPM 
Sol idity 


12/2.82 

100.8 

9.75 

10 % 

9 

2 

NASA 4418 (half rough) 
0 

34.7 

0.042 


Table III 
r/R 

Chord and 
Chord ft 

Twist Angle 

Twist Angle 
degree 

1 

2.82 

-3 

.9 

3.9 

-1.85 

.8 

5.0 

-0.62 

,7 

5.95 

0.62 

.6 

6.95 

1.85 

.5 

8.0 

3.0 

.4 

9.0 

4.28 

.3 

10.03 

5.5 

.2 

11.03 

6.69 

.1 

12.0 

8.0 



Fig. 1 - Force and velocity diagram for Glauert vortex theory. 
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Fig, 9 - Alternator power vs. wind speed; zero pitch angle 



Fig, 10 


Cp vs. tip speed ratio, MOD-1 (NASA 4415 airfoil] 



QUESTIONS AND ANSWERS 


D, Jeng 


From: Anonymous 

Q: Can you comment on actual power being higher than predicted at high wind speed? 

A: Comparison between experimental and theoretical values is valid only between wind 

speeds of 8 mph to 18 mph. When the wind turbine reaches the rated power of 100 kW , 
Because the blade geometry for the theoretical calculation was not changed ^ the 
actual power and predicted power cannot be compared for wind speeds which exceed 
the rated value of 18 mph (Fig, 9). 


From: Anonymous 

Q: Can any conclusion about aerodynamic correlation be drawn from cut in speed results 

which are largely dependent on mechanical friction and may vary from season to sea- 
son? 


A: 


To obtain 
losses by 
the text, 
machine is 
the actual 
the lowest 


alternator power^ the rotor power values were corrected for drive train 
an empirical equation applicable for the MOD-0 wind turbine and given in 
It is assumed that the power loss due to friction for this particular 
7,5 kW which is constant. By comparing the cut-in speed, we can compare 
rotor power (alternator power is zero) and the predicted rotor power at 
wind speed. 



